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Abstract 

The functions satisfying the mean value property for an n-dimensional cube are deter- 
^SJ ■ mined exphcitly. This problem is related to invariant theory for a finite reflection group, 

especially to a system of invariant differential equations. Solving this problem is reduced 
to showing that a certain set of invariant polynomials forms an invariant basis. After 
establishing a certain summation formula over Young diagrams, the latter problem is set- 
tled by considering a recursion formula involving Bernoulli numbers. 
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+^ 1 Introduction 

Let P be an ra-dimensional polytope in M". For = 0, . . . , n, let P{k) be the fc- dimensional 
skeleton of P. A continuous function / : — )■ M is said to be P{k) -harmonic if it satisfies 

S; f^^^ = T^(TV\l f{x + ry)dfik{y) (1) 

v/i^ , l^l'^JI Jp{k) 

in 

^ for any x G M" and r > 0, where is the fc-dimensional Euclidean measure on P{k) and 
T— I := l^k{P{k)) is the fc-dimensional Euclidean volume of P{k). This is an extension 

^ to a polytope of the classical notion of harmonic functions characterized by the mean value 
■ property for the (n — l)-dimensional sphere S*""^. Let 'Hp(fc) denote the set of all P(/c) -harmonic 
• I— I ! functions on R". A general result in Iwasaki [9j states that for any polytope P and any 
^ ' k = 0, . . . ,n, the set 'Hp(k) is a finite-dimensional linear space of polynomials. Note that 
' 'Hp(k) carries the structure of an M [5] -module, because equation ([1]) is stable under partial 
differentiations d = {di, . . . , (9„), where di := d/dxi is the z-th partial differential operator. 

It is an interesting problem to determine the space l-ip(k) explicitly when P is a regular 
convex polytope with center at the origin in R*^. This problem is already settled unless P is an 
n-dimensional cube. As for the cube case, however, although the vertex problem {k = 0) was 
solved by Flatto p] and Haeuslein [6] as early as 1970, the higher skeleton problem {k = 1, . . . ,n) 
has been open up to now. The aim of this article is to give a complete solution to this problem. 
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A characteristic feature of our work is that it provides a simultaneous resolution for all skeletons, 
which reveals a natural structure of this problem from the viewpoint of combinatorial analysis. 
Here we refer to Iwasaki [11] for a general review of the topic discussed in this article. 

Let C be an n-dimensional cube with center at the origin in the Euclidean space endowed 
with the standard orthonormal coordinates x = {xi, . . . , After a scale change and a rotation 
one may assume that the vertices of C are at (±1, . . . , ±1). The symmetry group of C is a 
finite reflection group of type which is the semi-direct product Wn := Sn ix {±1}" of the 
group {±1}" = {e = {ei, . . . ,5„) : Si = ±1} of n-tuple signs acting on by sign changes of 
Xi, . . . ,Xn, with the symmetric group Sn acting on by permuting xi, . . . ,Xn- The order of 
Wn is 2" ■ n! and the fundamental alternating polynomial of Wn is given by 

i<j 

The first main result of this article is then stated as follows. 

Theorem 1.1 Let C he an n- dimensional cube centered at the origin in M". For any k = 
0, . . . ,n, the linear space T-Lc^k) is of 2"- -nl- dimensions and as an M.[d]-module 'Hp{k) is generated 
by the fundamental alternating polynomial A{x) of the reflection group Wn- 

For an arbitrary polytope P Iwasaki [9j introduced an infinite sequence of homogeneous 
polynomials Tm\x) of degrees m = 1, 2, 3, . . . in terms of some combinatorial data about P{k) 
and characterized T-Lpik) as the solution space to the system of partial differential equations 

T^\d)f = (m = 1,2,3,...). 

From the way in which they are defined, the polynomials "* (x) are invariant under the sym- 
metry group G of P. This observation connects our problem to the theory of G-harmonic 
functions due to Steinberg [T?]. A C°°-function / : — )■ M is said to be G -harmonic if it satis- 
fies ip{d)f = for any G- invariant polynomial if{x) without constant term. Let T-Lg denote the 
set of all G-harmonic functions. There is always the inclusion Tic C l-ip[k)i and if {^'m' {x)}'^^i 
happens to generate the ring of G-invariant polynomials, then there occurs the coincidence 
T-Lg = ^p{k)- Steinberg [1^ made an explicit determination of T-Lg when G is a finite reflection 
group: dimT^G = 1^*1 and as an M[9]-module Hg is generated by the fundamental alternating 
polynomial of G. If P is a regular polytope, then G is a finite reflection group and we are done 
if we are able to show that the sequence {Tm\x)}^^i actually generates the G-invariant ring. 

For the /c-skeleton G{k) of the n-cube G the polynomials Tm\x) are constructed as follows. 
First recall that the m-th complete symmetric polynomial of j variables is defined by 

^^■nti,...,t,):= E (2) 

{mi,...,m,j)£'Pj (m) 

where Vj{m) is the set of all ordered partitions (mi, . . . , nij) of m by j nonnegative integers. 
Note that Hm^^\ti, . . . , tj, 0) = Hm{ti, . . . ,tj). For each k = 0, . . . ,n, we next define 

hl!^\x) := H^+^\xi + --- + Xn, X2 + ---+Xn, . . . , Xfc+i + ■ ■ ■ + Xn). (3) 

Note that when k = n the term Xk+i -|- • • • -|- a:„ is null and thus hm\x) = hm~^\x). 
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Figure 1: The cube in three dimensions with a flag V -< E -< F -< O 

For example, when n = 3 these polynomials are given by 

r H^\V-x) {k = 0), 

hl^\x) = \ H^^\V-x,E-x) (A; = l), 

[ Hi\v ■x,E-x,F-x) (A; = 2,3), 

where V = (1, 1, 1) is a vertex, E = (0, 1, 1) is the midpoint of an edge and F = (0, 0, 1) is the 
center of a face of the 3-cube C (see Figure [T]) and V ■ x stands for the inner product of V and 
X regarded as space vectors. Identify E and F with the edge and the face on which they lie. 
Similarly the origin O, i.e., the center of the cube C is identified with the unique 3-cell, i.e., the 
cube itself. Then we have a flag V -< E -< F -< O, where * -< ** indicates that * is a face of **. 
The simplex VEFO is a fundamental domain of the symmetry group W^. There is a bijection 
between the elements of and the flags of C. These pictures carry over in n dimensions. 
Finally Tm\x) is defined to be the PF„-symmetrization of hm\x), that is, the average: 

^i'n^):=^ E ^-(^^) ik = 0,...,n). (4) 
In other words Tm'' (x) is the average (x) over all flags of C. Note that 

since /irn^(^) ~ hrn ^(^) 

mentioned earher. It is immediate from definition that Tm ^ (x) is 
a homogeneous iy„-invariant of degree m. Recall that the degrees of Wn are 2,4, . . . ,2n, which 
are all even (see e.g. Humphreys [^). So the invariant polynomial vanishes identically 

for every m odd. The second main result of this article is then stated as follows. 
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Theorem 1.2 For any k = 0, . . . ,n, the polynomials t!^\x), Tf^\x)^ . . . , T2^^(x) form an in- 
variant basis of the reflection group Wn- 

For the proof of this we recaU that 62 (x), . . . , e2n{x) form an invariant basis of Wn, where 
e2m(a;) := ^ii ' ' ' (m = 1, . . . , n) 

l<n<---<im<n 

is the m-th elementary symmetric polynomial of x^, . . . , x^. Since r2^^ (x) is a homogeneous Wn- 
invariant of degree 2m, there exist a unique constant Cnjn and a unique weighted homogeneous 
polynomial Pn%{ti, . . . , tm-i) of degree 2m with tj being of weight 2i such that 



'2m 



(x) =ci%e2m{x) + Pi%{e2{x),...,e2m^2{x)) {m = l,...,n). (5) 



Here we employ the notation Cn,m and Pn,m. to emphasize the dependence upon n. Note that 
c^njn = oh,m^^ siucc T2^{x) = If we are able to show that the coefficient Cnju does 

not vanish for any m = 1, . . . , n, then we can invert equations (j5]) to express e2(a;), . . . , e2m{x) 
as polynomials of T2''\x), . . . ,T2m{x)- From this Theorem 11.21 follows immediately. So it is 

(k) 

important to develop a method to calculate Cn,m or at least to show that it does not vanish. 

(k) 

It turns out that the coefficients Cn,rn exhibit a beautiful combinatorial structure upon 
introducing the generating polynomials 

n , (fc) 

The third main result of this article is concerned with the structure of these polynomials. 
Theorem 1.3 The polynomials Gn,m{t) are tied to Gm{t) '■= Gm,m,{t) by a simple relation 

Gn,mit) = ii + ^T'^'GUt) (n > m > 1). (7) 
On the other hand the polynomials Gmif) admit a generating series representation 

Ti-ir-'Gut) (^Y = (8) 



m=l 



Equation ([8]) readily leads to a recursion formula for Gm{t) involving the Bernoulli numbers 
Bm- There are several conventions for defining Bernoulli numbers, but the most useful one in 
our context is through the Maclaurin series expansion 

00 „ 

^ — 1 _ £. I \ ^i' 1 \m-l ^rn ^2m 

e--l~ 2^ > (2m)! ' 

or equivalently through the formula 

z coth z = l + 2 Y.i-'^r-V^, bm := J^B„,. (9) 



m=l 



Multiplying formula (E]) by 2{tz coth 2 + 1), expanding the resulting equation into a power series 
of z"^, and comparing the m-th coefficients of both sides, we obtain the following. 
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Figure 2: An upper quadrilateral matrix 



Corollary 1.4 The polynomials G„i{t) satisfy a recursion formula 

m-l 

Gm{t) = b^{t + ir-' + 2t V + ir-'-'G,{t), G'i(t) = - + -. (10) 

^ — ^ z 

i=l 

A polynomial of degree m is said to be positive if its coefficients up to degree m are all 
positive. Note that the product of a positive polynomial of degree i and a positive polynomial 
of degree j is a positive polynomial of degree i + j. With definition we have 

C(2m) 1^1 
&- = ^ = ^E-^>0 (m = 1,2,3,...). 

i=i 

Therefore recursion formula flTU]) inductively implies that Gm{t) is a positive polynomial of 
degree m. Formula ([7]) then tells us that Gn,m{t) is a positive polynomial of degree n. Finally 
formula ([6]) concludes that the coefficient Cm,n is positive for any n > m > 1 and /c = 0, . . . , n. 
This establishes Theorem 11.21 The logical structure of our main results is this: 

Theorem 11.31 Corollary 11.41 > Theorem 11.21 > Theorem 11.11 

Thus the main body of this article is exclusively devoted to establishing Theorem 11.31 

The plan of this article is as follows. In Section |2] we represent the coefficient in terms 
of a sum over some matrices (see Proposition 12.51) . In Section [3] this representation is recast to a 
summation formula over some Young diagrams (see Proposition 13.41) . After these preliminaries. 
Theorem 1 1 . 3 1 and Corollary 11.41 are established in Section HI where some amplifications of these 
results and a summary on polyhedral harmonics for regular convex polytopes are also included. 



2 Matrix Representation 

(k) 

We derive a representation of the coefficient Cn,m as the sum of some quantities depending on 
a certain class of matrices. The main result of this section is given in Proposition 12.51 Various 
representations in this section involve those matrices as in Figure [H namely, A = (aij) with 
= for any i > j. Such a matrix is referred to as an upper quadrilateral matrix. Note 
that it becomes an upper triangular matrix if its vertical size is larger than or equal to its 
horizontal size. Throughout this article we use the following notation. For a matrix M = {mij) 
of nonnegative integers, whether upper quadrilateral or not, or even for a row or column vector. 

Ml := J]m,,!. 

hi 
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For a row vector v = {vi, . . . , f„) of nonnegative integers we put := x^^ ■ ■ ■ x^. Moreover, 



fc+i 

e := ('1,^ r), 



1 := 



/1\ 



Lemma 2.1 The polynomial in is expressed as 



X 



where the sum is taken over all {k + 1) x n quadrilateral matrices A of nonnegative integers 
whose entries sum up to m. 

Proof. In view of definitions ([2]) and ([3]), the multi-nomial theorem yields 



k+l 



mi-\ [-"1^+1=771 i=l 

k+l 



E n E rn^U', 

miH |-mfc+i=m, j=l XaaA hai„=mi •'■ •'■i— * j=i 



J2aij=m 

E 



H h a»n)! 

ni<, « 



i<j "'V 



n 

i<j 

n 

n 



Emin{ j,k4-l\ 



where atj is defined for l<i<j<n, i<k + l. Putting aij = for k + l>i>j>l makes 
A = [oij) an upper quadrilateral matrix. It is obvious that the entries of A sum up to m. □ 

Consider the {±l}"-symmetrization of hm\x), that is, the average: 



4 Yl hl^\^iXi,---,enXn), e =(£!,...,£„) e {±1}". 



(12) 



£G{±1}" 



Lemma 2.2 The polynomial in (ITlj) is expressed as 

^^'^(-) = E^-'^ (13) 

A 

where the sum is taken over all {k + 1) x n quadrilateral matrices A of nonnegative integers 
whose entries sum up to m and moreover whose column-sums are all even. 
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Proof. Substituting formula ( 1TT|1 into definition f[T2l) yields 

££{±1}" A A \ee{±l}" J 

where the matrix A ranges in the same manner as in formula flTT]) . Put eA = (z/i, . . . , i/„), 
where uj is the j-th column-sum of A. Observe that 



E,-^ I 2" (uj is even for any j = 1, . . . , n), 

£€{±1}" ee{±l}" I J = 1, • • • , 



So the sum in (|T4l) can be restricted to those A's whose column-sums are all even. □ 

For any matrix with even column-sums, its entries must sum up to an even number, so 
that formula f[T^ implies that gm\x) vanishes identically for every m odd. Thus from now 
on m is replaced by 2m with m being a positive integer. This allows us to put eA = 2i'{A) 
with z/(A) = (z/i(A), ...,!/„( A)) G Vn{fn), that is, i^iA) is an ordered n-partition of m. The 
polynomial Tg^l^) formula (jl]) is the Sn-symmetrization of gimi^)^ that is. 



n! 



^ XI ^'2m(a^a{l).---.2;^(n))- (15) 



Putting formula (|T3l) with m replaced by 2m into formula (fT5|) yields 

>P (^1)! 2^1 (M) 2i.n(A) ..^x 

where J^n}m is the set of all (fc + 1) x n upper quadrilateral matrices of nonnegative integers 
whose entries sum up to 2m and moreover whose column-sums are all even. 
Let denote a primitive j-th root of unity. Since 



e2j(C2m5 ^2mi • ■ • i ^Imi 0, • • • , 0) 



(j = l,...,m-l), 

(-ir-i (j = m), 

substituting x = (C2m, Clm^ • • • 5 C2m5 0, . . . , 0) iuto equation (jS]) yields 

= {-ir-'4'2{C2m, CL, . . . , C2^, 0, . . . , 0). (17) 

For each partition z/ = (z^i, . . . , z/„) G Vn{rn), we define 

w„,^(j^) = n„(z/i, . . . , z/„) := Yl C^^^^^+-+'^(")^", (18) 

where Sn,m{^) '■= {o' E Sn '■ for i = 1, . . . , n, if z/j > 1 then a{i) G {1, . . . , m} }. Since 

(a G 5„,„(z.)), 



. . . „2z/„ 
^<7(1) ^a{n) 



(o- ^ ^^^^(z/)). 



at a; = ((2™, Clm^ • • • 5 (2^^ 0, . . . , 0), formulas ([IS]) and ((IT]) yield 



E (19) 



Let denote the number of positive entries in an ordered partition v = [ui, . . . , z/„) G 
Vnijn). Note that £(z/) G {1, . . . , m}, because ^(z/) < z/i + ■ ■ ■ + z/„ = m. 

Lemma 2.3 The function Un,m{v) is symmetric, that is, invariant under any permutation of 
1^1, ... , z/„. For any v = {ui, . . . , G Vn{rn) such that I'm+i = ■ ■ ■ = = 0, we have 

[ji _ £(z/))! 

Wn,m('^) = 7 7^r^Mm,m(z/i,...,Z/m). (20) 

Proof. For an element r G put z/"^ := {vrii), ■ ■ ■ , ^rin))- Then it is easy to see that Sn,m{^^) = 
Sn,m{j^) ■ T- Using this we show that Mn,m('^'') = Un,m{y)- Indeed, 



as desired. This proves that Un,m{^) is a symmetric function oi v = (i^i, . . . , z/„). 

We proceed to the second assertion. Suppose that v is of the form v = (i^i, . . . , z/,,, 0, . . . , 0) 
with r := C{u) < m. Then Sn,m{y) = {o" G S'„ : cr({l, . . . , r}) C {1, . . . , m}}. We think of Sm 
as a subgroup of by setting Sm '■= {cr ^ Sn '■ crii) = i for i = m + 1, . . . , n}. Define a map 

a{i) (2 = l,...,r), 

5'„,m(z^) ^ Srn, (T ^ T by r(i) := { p(i) (z = r + 1, . . . , m), (21) 

i (z = m + 1, . . . , n), 

where p is the unique bijection p : {r + 1, . . . , m} — )■ {1, . . . , m} \ a({l, . . . , r}) which is "order- 
equivalent" to the injection cr\^r+i,...,m} in the sense that p{i) < p{j) if and only if a{i) < cr(j) 
for every i,j G {r + 1, . . . , m}. We claim that the map -one. Indeed, given any 

element r G S*™,, the fiber over r has a one-to-one correspondence with the set of data {S, q): 

• a subset S of cardinality m — r of T := {1, . . . , \ r({l, . . . , r}), 

• a bijection q : {m + 1, . . . , n} — )■ T \ 5. 

It is clear from definition fl^Tl) that given a data (S*, q) there exists a unique element a G 5'„^m(i^) 
such that cr({r + 1, . . . ,m}) = S and cT|{m+i,...,n} = ^Z- Since #T = n — r, there are (^l^!,) 
choices of S", for each of which there are (n — m)! choices of q. Thus the fiber has a total of 
(r;) ■ - ^)! = ^ elements. Since = ^rji)u,+-+rMu^ ^ have 

(n-r)! ^ (^-0' ^ ^ ^ 

(m — rj! \m — r)\ 

where Sm = Sm,m{^i, ■ ■ ■ , ^m) is used in the second equality. □ 
Formula ( 12U]) reduces the calculation of m„ ^ to that of Um,m, which we now carry out. 
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Lemma 2.4 For any partition v = (z/i, . . . , Vm) € Vml^m), 

u.m{v) ■■= UrnA^) = m (-1)^^-1 - 1)! (m - £(z/))!. (22) 

Proof. When n = m the function Un,m{i^) in definition f|T8|) becomes simpler because Sm,m{i^) = 
Sm for every u G Vm{m). The proof is by induction on £(z/). When = 1, we may assume 
that u is of the form u = (m, 0, . . . , 0) by the symmetry of ^^(z^). Then definition ([T5|) reads 

= E C^'^" = 5^ 1 = m!, 

which verifies formula ( 12 2 p for l{v) = 1. Let 1 < r < m and assume that formula (1221) is true 
for every partition u G Vm{m) with £(z/) = r. Consider the case £(z/) = r + 1. By the symmetry 
of Um{i') we may assume that u is of the form z/ = (z/i, . . . , i^r+i, 0, . . . , 0) with i/i, . . . , z/,.+i > 1 
and z/i + ■ ■ ■ + z/,.4.1 = m. Note that 1 < Vr+i < Formula (ITSl) now reads 



O-gSm (pi,...,Pr+l) 

where {pi, . . . ,Pr+i) ranges over all permutations of distinct r + 1 numbers in {1, ... , m}. Thus, 

(pi,...,Pr) pr+l£{l,...,m}\{pi,...,pr} 



'-Pr + l'^r+1 

5m 



(pi,...,Pr) V« = l J' = l 

Since 1 < z/^+i < we have X]/=i Cm"^^^ = and hence 



z — / z — / m — r 

j=i (pi,...,pr) i=i 

where z/^-^^ = (z/p\ . . . , I'm) := (z/i + SijUr+i, . . . , z/i + drjVr+i-, 0, . . . , 0) with 5ij Kronecker's 
symbol. Note that for each j = 1, . . . , r, we have z/'^-'^ G Vm{^) with r(z/(-^^) = r, so that the 
induction hypothesis yields Umit^^-'^) = rn (—1)''^^ [r — \)\ [m — r)\ for j = 1, . . . , r. Therefore, 

Um.{^) = f ("l)*" ""^ {f^ ~ 1)! (m — r)! = m (— 1)^ r! (m — r — 1)!, 

m — r 

which means that formula (1221) is true for £(z/) = r + L The induction is complete. □ 
A column of a matrix is said to be nontrivial if it has at least one nonzero entry. 

Proposition 2.5 Let i{A) denote the number of nontrivial columns in A. Then, 

Proof. First, Lemmas 12.31 and 12.41 are put together to yield the formula 

UnA^) = m (-l)^('^)-^ (£(z/) -l)\{n- £(z/))!, (24) 

for any partition u = (z/i, . . . , z/„) G Vnijn). Indeed, by the symmetry of Un,m{,^) "^^ assume 
Vm+i = ■ ■ ■ = i^n = 0. Thus using formula (1221) in formula ( 12U]) gives formula (I21|) . Next, putting 
formula 021]) with z/ = z/(A) into <^ yields formula ([2S]), since i{A) = ^(z/(A)). □ 
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3 Young Diagram Representation 



We rewrite formula (123 p in Proposition 12.51 as a sum over some Young diagrams. After several 
preliminary discussions, the main result of this section is stated in Proposition 13.41 For each 
V = [ui, . . . , z/„) G Z>Q, let Ain \i') be the set of all [k + 1) x n upper quadrilateral matrices 
the i-th column of which sums up to z/j for i = 1, . . . , Motivated by expression fl23l) . put 

-fM:= E (25) 



Lemma 3.1 For any v = (z/i, . . . , z/„) G Zi'>q, we have 

n-=i(^i + ---+^.+j)-n;=i^.! 

Proof. The proof is by induction on k. For A; = there is nothing to prove. Suppose that 
formula (126|) is true for — 1. We write v = v to emphasize that i/ is a row vector. Put 

?/'^''^ (ai, . . . , 4+i) := -^^^p for A = : with = (^7"^^, a^i, . . . , aj„). 

Here we also write dk+i = a = (0, . . . , 0, ak+i, ■ ■ ■ , an) to simplify the notation. Observe that 
il)n\ai, . . . ,ak,a) = ipn\a) ■ ipn~^\ai, . . . , Sk). Using this we have 



^ ijit\ai,...,ak,a) 

diA \-ak+a=u 



a<£' 



ffliH \-ak=u—a 

where a < U means that u — a E Z>q. Put /i^ := Ui + ■ ■ ■ + Uj, p,j := Uj^i + ■ ■ ■ + z/„ and 
b := afc_|_i + ■ ■ ■ + a„. Since = for j = 1, . . . , fc, the induction hypothesis yields 

1 + 



Substituting this into the previous formula and after some manipulations we have 

1 



n,t?(/^.+j)-n;=i^.! 



E&!(^n + A:-i-6)! E n 



6=0 



E&! (/in + A;-l-6)! 



b=0 



d<u j=k+l 
o.k+1 + ■ ■ ■ + a„ = b 



b 



fin + k - 1 - b 
Hk + k-l 



/ifc! (/ifc + A: - 1)! 

n,t;(/^.+j)-n;=i^.!t^ 

j^fc! {jJik + k- 1)! f^in + k 

iiYJM + j) ■ n;=i U + kj n .tii/i, + j) ■ n;=i ^ 
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c, 



where the following general formulas are used to obtain the second and fourth equalities. 

0<b,<a,i=l^^ ^ ^ i=b ^ ^ ^ ' 

bi + ---+bk = b 

Therefore formula (126|) remains true for k and the induction is complete. □ 
Lemma 3.2 Formula ( \23^ in Proposition \2.5\ is rewritten as 

<^= ^'^^"'"'!^^^- E {-lY'''-\^{^)-l)Kn-tW))\€\-\ (27) 

where 

vi^Kj') ■= • (28) 

n-=i(2^i + ---+2z.,+j)-n;=i(2^,)! 

Proof. Since there exists a direct sum decomposition 

and £{A) = i{v) for A e Mii'\2u) with u E V„{m), formulas ([23]) and definition ([25]) lead to 
= E (-1)'""-'W-) -!)!("-%))! E ^ 

= ^"'^7'" E (-ir^^^-^ (^(-)-l)!(n -£(.))! (2.). 

Use formula fl26|l with z/ replaced by 2z/ and factor the term (2z/i + ■ ■ ■ + 2z/„ + fc)! = (2m + k)\, 
which is constant for u G Vn{m), out of the summation. Then we obtain formula (12 7p . □ 

Formula (I27|) comes up as a sum over ordered partitions. The next task is to recast it to a 
sum over unordered partitions, that is, over Young diagrams. Let y,- be the set of all weakly 
decreasing sequences A = (Ai > ■ ■ ■ > A-,) of nonnegative integers. The sum |A| := Ai + • — \- Xj 
is called the weight of A. Note that A represents an unordered partition of |A| by j nonnegative 
integers. The number of positive entries in A, denoted ^(A), is called the length of A. An element 
A G defines a Young diagram of weight |A| and of length i{X) < j in the usual manner (see 
e.g. Macdonald [IS])- An element A G is also written A = (0^'' l*"^ 2'"^ ■ ■ ■ ) when the number 
i occurs exactly times in A for each i = 0, 1, 2, ... , where the term z*"' may be omitted if 
rj = 0. Note that |A| = ri + 2r2 + Sr^ H — ■ , £{X) = ri + r2 + r^ + ■ ■ ■ , and j = To + ^(A). Given 
an element A G yj let Vj{X) denote the fiber over A of the order- forgetful mapping 

^>o ^ ^ = (^i> . . . , z/,) ^ A = (0^° ■■■), 

where i occurs Sj times in u for each i = 0,1,2, ... . Motivated by expression (1281) . consider 

Wk(u) := V ^ r (29) 

.,tt,)n,U2.i + ---+2.,+j)-n,ti(2..)! 

for fi E where the denominator of the summand differs from that of formula (1281) by the 
factor rij=i(2^'j)! in place of ]n!j=i(2i^j)!- This function is evaluated in the following manner. 
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Lemma 3.3 For each /j, = {0^^° V^^ 2^*2 ■■■ ) g yt, we have 

Proof. The proof is by induction on k. For k = 1 formula ( 13 Up holds trivially. Suppose that 
k > 2 and formula f l5U]) is true for fc— 1. Let = (j''^ | j G J) G 3^^ with J := {j G Z>o : Sj 7^ 0}, 
where ^jgjjsj = and ^jgjSj = /c. For each i G J, put /i*^') := \j G J), where 5ij 

is Kronecker's symbol. Since /i'-*-* G 3^fc_i, the induction hypothesis implies that for each i G J, 

Observing that there exists a direct sum decomposition 

and noticing that 2z/i + ■ ■ ■ + 21/^ + k = 2\fi\ + k, we have 

° £ (2|H + . (2,)! nS(2.!" + ■ ■ ■ + 2.f + nS(2^")! 

= g (2|H+l).(2.)! "'=-''^'"' 

1 £i6j-^i(2i + l) 

(2|f.| + *:) n,«^j!n,«((2i + l)!)'' 

This shows that formula (130|) remains true for /c and the induction is complete. □ 
Let 3^n(^) be the set of all unordered n-partitions of m and put y{m) := 3^m(m). 

Proposition 3.4 The generating polynomial Gn,mit) in' definition (0) is expressed as 

^^^^> xmm) n.---r^. 

for any n > m > 1, where A = (0^'' ■ ■ ■ m'"™) and Tj is defined by 

^ (2j + l)! t+l u > > ; \ ) 

In particular the rational function {t + l)^"'Gn,m{t) is independent of n. 
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Proof. For each A G 3^n("^) let 3^fc(A) be the set of all Young subdiagrams yU G 3^^ of A. For 
each fi G 3^a:(A) let Vn\\, ^) be the set of all z/ = (z/i, . . . , G "PnlA) such that the cut-off to 
the first k components u' := (z/i, . . . , i/^) belongs to Vkifi)- Let A = (C" 2*^2 ■ ■ ■ ) G ^^(m) 
and /i = (0'*" •■■) G 3^fc(A). Note that rj = sj = for j > m. Taking /i away from 

A induces a skew-diagram A//i = (0^°"*° xn-^i 2''2-«2 . . . ^ with rj — Sj > for j = 0, 1, 2, . . . . 
Since #P„-,(A/m) = (n - A;)!/n,>o(^. - ^.)! and Utk+ii^^iV- = U,>oiiW-p-'^ , we have 

y v(^\u) = y ^ 

l<j<fc l<i<A: j>0 j>0 



n,>o^.!yo W V(2j + !)!>' V(2j) 
Since there is a direct sum decomposition P„(A) = U„g-w Vn\\,ij), we have 



T— — y #)(z/) = v— — y y v^'^Hu) 



so + ■^l + ■ ■ ■ — 



n."Lo^.!S'"£oPoV^.A(2j + l)!; V(2i)! 

1 / 1 t 



nr=o-.!|iV(2j + l)! (2j)! 

Jl + 1)' 
ro! 



m 



n^'^ (34) 



where n — k = J2f=oi^j ~ ^j) ^ — Sj=o '^i used in the third and final equalities. 
On the other hand, in view of Vn{m) = Uxeynim) '^nW, formula (I27|) yields 
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Thus, taking i{X) = ri + ■ ■ ■ + and n — i{X) = ro into account, we have 

(_l)m-i'^nAt) ^ (-iTl^A n\ cgL 1^ (by m) 

^ ^ + {t + ir ^ {n-k)\{2m + k)\ 



k=0 



(t + 1)" ^ ' ' ' ' ' ^ ^ - A;)! 



m 

Ae:i^„(m) ™ 



where the sum may be taken over 3^(m), because when n > m any A = (Ai > ■ ■ ■ > A„) G 3^n('^) 
is of length at most m, that is, Xj = for any j > m, and hence A can be identified with 
A' := (Ai > ■ ■ ■ > Am) G 3^(m). This proves formula fl52]) . As the right-hand side of formula 
( 132|) depends only on m, the rational function (t + l)''"'Gn,m(t) is independent of n. □ 



4 Generating Functions and Bernoulli Numbers 

We are now in a position to establish Theorem 11.31 and Corollary 11.41 

Proofs of Theorem \1.3\ and Corollary I.4 Formula ([7]) is an immediate consequence of the last 
assertion in Proposition 13.41 that {t + l)~"'Gn,m{t) is independent of n. The proof of formula 
dH]) is based on the following general fact on generating series: if we put 

(3m ■.= m J2 (-l)-+-+--ii!l±:-:±l!IL^ a? ■ ■ ■ a- (m = 1, 2, 3, . . . ), 

ri! ■ ■ • Tm' 

where A = (0^° 1*^^ ■ ■ ■m''"'), then there exists a formal power series expansion 

log(l + f:a..-)=f;^.- (35) 

\ m=l / m=l 

We apply this formula to the situation of Proposition 13. 4[ where ctj = Tj in formula ( 133|1 and 

/3„ = (-l)'"-i(t + l)-'"G'm(t) (36) 
in formula fl5^ with n = m. We now find 



^ (2m)! t + l^(2m + l)! 

1 ^ ' m=l ^ ' 



771=1 



^ ^ 1 ^ 

T ^ f2m)! ^ t+1 5^ 



t + 1 ^ (2m)! t+l^(2m + l)! 

■m=0 ^ ^ 7Tt=0 ^ ' 

sinhzA . 
tcosh;2 + . (37) 



t + 1 

Substitute this into formula (l35i) and apply the differential operator f ^ to the resulting equa- 
tion. Then after some calculations we get formula (E]) and thus establish Theorem II. 3 [ Corollary 
11.41 then follows easily from this theorem in the manner mentioned in the Introduction. □ 
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We present some amplifications of Tlieorem 11.31 and Corollary 11.41 For the extremal cases 

(k) 

of /c = 0, n — 1, n, the coefficients Cn,m can be written explicitly in terms of Bernoulli numbers. 



Lemma 4.1 For k = 0, n — 1, n, the coefficients Cnjn are directly tied to bm by 

c'^l={2m)\{2'"^-l)b^ c(::-^^ = cW = (n>m>l). (38) 

Proof. Substitute t = into definition (Q to get G'n,m(0) = i^n+2m)\ '^""m- Similarly put t = 
in formulas ([7]) and ffTOj) to have G„,m(0) = Gm{0) = b^- These together lead to the assertion 
for ci"m in formula ( 138|) . The assertion for ci"m^'' then follows from the identity Cn"m^'' = ci"m 
mentioned in the Introduction. To prove the assertion for Cn}n in formula f l38p we consider the 
generating polynomial Gn,m{'t) '■= t"'Gn,m{^/t) instead of Gn,m{t). After the change t (-)■ 1/t 
and multiplication by t", formula (jH]) gives G„,m(0) = Cn°L/(2m)!. On the other hand, formula 
([7]) yields Gn,m{'t) = + ^Y'^Gmit), where Gm{t) '■= Gm,m{t), while formula ([HD gives 



°° / 2 

5^(-l)-^G.(t)(^ 

m=l ^ 



t(2;cothz — 1) + 
2{zcoth z + t) ' 



which upon putting t = reduces to the equahty 



' ^ ^ 2cothz 2 

m=l 

Comparing it with the Maclaurin expansion tanh^; = 2 ^^^-^(— 1)™^"'^(2^™ — l)bmz'^"^~^, we 
find G^(0) = (22™ - 1) bm. Thus = (2m)! ^^.^(O) = (2m)! G„(0) = (2m)! (2^™ - 1) □ 

The first formula in ( l38|) is already found in To deal with the intermediate coefficients 
Cn)n for /c = 1, . . . , n — 2, another modification of the generating polynomials Gn,m{t) is helpful. 

^)-E („_^,;;:^,), '^(i-')- (39) 



Lemma 4.2 For n > m > 1, t/ie polynomials Fn^rn{t) depend only on m, being independent of 
n. They satisfy the differential- difference equation 

2F„,™(t) + -F;„(t) + il^F^,^_,(t) = (n > m > 2). (40) 
m ' m — 1 

All the Fn^rn{t) Can be determined inductively by solving equation with initial conditions 
Fn,m(0) = (2^™ - 1)6^, Fn,i{t) = ^-^ {n>m>l). (41) 
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Proof. Put Fm(t) := Fm,m(t)- It readily follows from relation ([7]) and definition fl39l) that 
Fn,m{t) = Fm{t) for every n>m. The substitution 1 1-)- ^ induces the changes 

°° ■ u 

H-> (-1)— 1 + ^ r™;^^"^ ^ (1 - t) coshz + t 

m=l 

in formulas f l36p and (I57|) respectively. With these changes formula fl35p reads 

log ((1 - t) cosh^ + t = izl)!liF^(t)^2™. (42) 

^ m=l 

Denote the both sides of this equation by $ = $(2;, t). A direct check using the left-hand side 
of equation tells us that $ satisfies the partial differential equation 

+ {t-(l-t)V} — = (l-t)^l (43) 

Next we look at this equation by means of the right-hand side of formula fH2|) . For each m > 2 
the coefficient of z^™ in equation fj^3l) being zero gives the differential-difference equation 



2FUt) + -F^t) + ^-^-^Fi^,{t) = (m > 2), 

which can be expressed as equation fjlO]) . because -Fm(t) = Fn^mit) and Fm-i(t) = 
by the first assertion of the lemma. The first condition in (14T|) is derived from formulas (l39i) 
and fl38l) as -F„_m(0) = ci°L/(2m)! = (2^"^ — l)bm, while the second condition follows from 
Fn,i{t) = Fi^i{t) and the direct calculation of Fi^i{t), which is easy. □ 

Differential-difference equation ( HUj) can be used to derive a recursion formula for Cnjn as 
well as to explicitly determine c^)n for k near or n in terms of Bernoulli numbers. 

Proposition 4.3 For k = 0, 1, the coefficients Sn}m are given by the first formula in f[3'5|) and 
c(^|„ = (2m + 1)! |(22- - 1) 6„ - _ 1) h^+i^ (n > m > 1). (44) 

For k = n ~ 2, n — 1, n, the coefficients Cn}m take a common value which is given by 

'^n.m ^ ~ '^n.m ^ ~ '^n.m ~ j " {u > m > 2) . (45) 

nl 



Moreover for 1 < k < n ~ 2 and 2 < m < n there exists a recursion formula 

c^-^^ = "'''^"^ hm + k-l) -ik + 1) ct\lJ . (46) 

nim — 11 ' 'J 



Proof. Write the left-hand side of equation ( HQ]) as X]fc=o ^"^"^ t'^(l — t)"" Since the polynomials 
t^{l — 1)"~'^, k = 0, . . . ,n, form a basis of the linear space of polynomials in t of degree at most 
n, it follows from equation ( HOj) that 7!^^ = for every < A; < n. For A; = we find 



7: 



uL = H T < 7:^ — ^-TTT - — — > = (n>m>2) 



(2m)! m-l|(2m-l)! (2m - 2)! 
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Figure 3: Approximations of the sphere by geodesic domes 



where Cn,™. and Cn\,^_i are aheady known as in the first formula of fl38|) . Thus c„_i is 
also known from this equation. Replacing {n — l,m — 1) with {n,m) we get formula (jS]). 
On the other hand, for k = n — 1, n, some calculations show that 7!"™^'* = and 7^"™ = 
lead to c^n,m'^ = c"n,rn^ and ci"m^'* = Cn"m respectively, where the latter is already pointed out 
in the Introduction and Lemma 14.11 Thus formula f H5]l follows from the second formula of 
f l35]l . Finally some more calculations of 7m,n for general k imply that equation '-ynjn = with 
1 < k < n — 2 and 2 < m < n is equivalent to recursion formula ( H6|) . □ 

Since c^Jn is already known for the fc's at both ends of the interval < A; < n as in 
formulas ( 138|) . f H4|) and ( H5|) . the recursion formula (146|) can be used to inductively determine 
all coefficients Cn,m, where there are three directions in which induction works productively. 

/ N (A:) (fc-l) (A:) (fc+1) /in (fc-1) (fc) (fc) (fc+1) 

l^aj Cn,m ^ Cn,m , C„_x,m-l5 '^n-l,m-l5 l^j Cn,m <r- Cn,m, C^-l^m-l^ ^n~l,m-li 

(c) c'^.^lln^i <r- c^n)m-, ci^m^\ c^'^^^m.i (with (m - 1, ^ - 1) replaced by {m,n)). 
For example formula f H6|l with k = n — 2 is used in direction (b) to derive 

>^3) ^ ^ 2m)! b^-4m- (n + 2m- 3)! (n > 3, n > m > 2) 

from formula fHSjl . Similarly formula f H6|l with = 1 can be applied in direction (c) to deduce 
a closed expression for Cm,n from formulas ( 138|) and f H4|) . and so on. 

At the end we return to the starting point of this article, that is, to polyhedral harmonics. 
With Theorem 11.11 for the cube case, the determination of polyhedral harmonic functions for 
all skeletons of all regular convex polytopes has been completed. As a summary we have: 

Theorem 4.4 Let P be any n- dimensional regular convex polytope with center at the origin in 
and G the symmetry group of P. Then for any k = 0, . . . ,n, the linear space l-ip^k) is of 
\G\- dimensions, where \G\ denotes the order of G, and as an M\d]-module 'Hp(k) is generated 
by the fundamental alternating polynomial Aq of the reflection group G. 

For the classification of regular convex polytopes we refer to Coxeter ^Lj. Theorem 14.41 is 
proved in article [lOj for the n-dimensional regular simplex and in article [12] for the exceptional 
regular polytopes, that is, for the dodecahedron and icosahedron in 3-dimensions and for the 
24-cell, 120-cell and 600-cell in 4-dimensions. For the n-dimensional cross polytope, namely, the 
analogue in n-dimensions of the octahedron, there is no detailed written proof in the literature, 
but a proof quite similar to the regular ra-simplex case is feasible. This is because each face 
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of an n-dimensional cross polytope is an {n — l)-dimensional regular simplex. Finally the n- 
dimensional cube case has been treated in this article (Theorem ll.ip . in which case the proof 
is quite different from those in the other cases. Here we should also mention the important 
studies [21 El in El El E] etc. of earher times, which contain partial answers to our questions, 
referring to the survey [TT] for a more extensive literature. 

Apart from the regular figures for which symmetry plays a dominant role, polyhedral har- 
monics is largely open, for example, for such figures as geodesic domes in Figure El 
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